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We study analytic structure of the Green's function (GF) for the exactly solvable Fano- Anderson 
model. We analyze the GF poles, branch points and Riemann surface, and show how the Fermi's 
Golden Rule, valid in perturbative regime for not to large time, appears in this context. The knowl- 
edge of analytic structure of the GF in frequency representation opens opportunities for obtaining 
formulas for the GF in time representation alternative to the standard one using the spectral density. 
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The Fano- Anderson model iJj], which presents discrete 
level coupled to continuum is probably one of the simplest 
and best known in quantum mechanics. It is exactly solv- 
able, and the solution of the model in terms of Green's 
functions (GF) is presented in The spectral line in- 
tensity calculated in the paper by Fano are conveniently 
represented using the Green's function spectral density 
• When the problem is treated within the quantum me- 
chanics proper, one is typically interested in the tunnel- 
ing of the particle, initially localized at the discrete level 
into contimmm (see e.g. Cohen- Tannoudji et. al. [^). 
The non-tunneling amplitude is standardly expresses also 
through the Green's function spectral density, so prag- 
matically speaking this is what we need. However, in 
this paper we would like to see the Green's function in a 
broader context, and study its analytical structure as a 
function of frequency. We'll see that Green's function (in 
frequency representation) is a multi valued function, and 
on simple examples study its branch points, poles and 
Riemann surface. In fact, we'll see that this study can 
be of practical value also, because it opens new oppor- 
tunities to connect the frequency representation of the 
Green's function (used in calculation of the spectral line 
intensity) with the time representation, necessary for the 
quantum mechanics proper. 

To have everything at hand let us present a pedagog- 
ical derivation of the non-tunneling amplitude. Our sys- 
tem consists of the continuum band, the states bearing 
index fc, and the discrete state d, having energy e. The 
Hamiltonian of the problem is 



H = Y,^k \k) (fc| + e \d) {d\ + {Vk \k) {d\ + h.c), (1) 

fe 



where |fc) is a band state and \d) is the state localized 
at site d; h.c. stands for the Hermitian conjugate. The 
wave-function can be presented as 



m - git) \d) 



k 



b{k,t)\k), 



(2) 



with the initial conditions g{0) = 1, b{k,0) = 0. Notice 
that the non-tunneling amplitude is just the appropriate 
GF in time representation. Schroedinger Equation for 
the model considered takes the form 
dg{t) 



dt 



Making Fourier transformation (Im a; > 0) 



(3) 



(4) 



we obtain 



-i + ujg{uj) = eg{u;) + YVkb{k,uj) 

k 

ijjb{k,uj) = LUkb{k,uj) + Vkg{uj). 



(5) 



For the amplitude to find electron at the discrete level, 
straightforward algebra gives 
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g{u})e-"^'duj 



where 



and 



uj — e — S(cij) ■ 
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UJ ~UJk 



(6) 



(7) 



(8) 



The integration in Eq. ([6]) is along any infinite straight 
line parallel to real axis in the upper half plane of the 
complex Lu plane. Notice that (?(a;) is the GF in frequency 
representation. The quantity T,(uj) is self-energy (or mass 
operator). 

For tunneling into continuum, the sum in Eq. ^ 
should be considered as an integral, and Eq. ([5]) takes 
the form 



S(u.) 



A(g) 



dE, 



where 



A(i?)=^|Ffcp<5(i?-a;fe), 



(9) 
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where and the limit of integration are the band bottom 
Eb and the top of the band Et . We would like to calculate 
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where 



FIG. 1: Contour used to evaluate integral (|6}. Radius of the 
arc goes to infinity. 



l/r = 27rA(e) 



(17) 



However, even in the perturbative regime, the FGR has 
a limited time-domain of applicability [3] . For large t the 
survival probability is determined by the contribution to 
the integral coming from the end points. This con- 
tribution can be evaluated even without assuming that 
the coupling is perturbative. Let A{E) ^ [E — Eb)^ 
/3 > 0) near the band bottom (the contribution from the 
other end point is similar). Then for large t 



(18) 



The similar contribution comes from the top of the band. 

If there can exist poles of the propagator ([7|) , we should 
add the residues to the integral (HH). Thus we obtain 



integral ([6|) closing the integration contour by a semi- 
circle of an infinite radius in the lower half-plane. Thus 
we need to continue analytically the function giuj) which 
was defined initially in the upper half plane (excluding 
real axis) to the whole complex plane. We can do it 
quite simply, by considering Eqs. (O and (O as defining 
propagator in the whole complex plane, save an interval 
of real axis between the points Ei, and Et , where Eq. ([9]) 
is undetermined. (Propagator analytically continued in 
such a way we'll call the standard propagator.) Thus the 
integral is determined by the integral of the sides of the 
branch cut between the points E^ and Et- The real part 
of the self-energy S' is continuous across the cut, and the 
imaginary part S" changes sign 



T."{E + is) = £"(£;- 



TrA{E) s +0. (11) 



So the integral along the branch cut is 

"^^ A{E)e-'^'dE 
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E, [E-e-^'{E)Y +^^A'^{E) 



Thus we have 

g{t) ^ icutit), 

and the survival probability p{t) is 
P{t)^\9{t)\\ 



(12) 



(13) 



(14) 



In the perturbative regime |E'(e)|, |S"(e)| <C e — 
Eb, Et—e the main contribution to the integral (fT^ comes 
from the region E ^ e. Hence the integral can be pre- 
sented as 



Icut — 



A{e)e-'^*dE 



(£;-e-S'(e))2+7r2A2(e) 



(15) 



and easily calculated to give the well known Fermi's 
golden rule (FGR) 



Pit) 



(16) 



git) ^lcutit) + 



(19) 



where the index j enumerates all the real poles Ej of the 
integrand, and 



-iEit 



Hi 
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dE \E=Ei 



(20) 



Notice, that the poles correspond to the energies of bound 
states which can possibly occur for E < Ef, or E > Et, 
and which are given by the Equation 



E, 



Set 

k ■' 



(21) 



If we take into account that normalized bound states are 

M> + E.^|fc> 



\E,) 



1 



1/2 ' 



(22) 



then the residue can be easily interpreted as the ampli- 
tude of the bound state in the initial state times the 
evolution operator of the bound state times the ampli- 
tude of the state \d) in the bound state 



R, = (d \E,) {E, \d) , 



iEjt 



(23) 



If the propagator has one real pole at Ei, from Eq. (fT9|) 
we see that the survival probability p{t) — > when 
t oo. If there are several poles, this equation gives 
Rabi oscillations. 

Notice that Eq. is just the well known result 



where 



lit) 



A{uj)e 



-Im [g{E + is)] 



(24) 



(25) 
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is the spectral density function. The first term in Eq. 
(fT9)) is the contribution from the continuous spectrum, 
and the second term is the contribution from the discrete 
states. 

An indication that there is more in the GF than we 
have so far discussed comes from the following fact: we 
could have obtained the FGR in perturbative regime di- 
rectly from Eq. ([6]), changing exact Green function ([7]) 
to an approximate one, which may be called the FGR 
propagator 



1 



S'(e) +i7rA(e)' 



(26) 



Thus approximated, propagator has a simple pole uj = 
e + S(e), and the residue gives Eq. (fT6|) . Thus, the 
propagator which is used to obtain Eq. (jlSp and the FGR 
propagator have totally different singularities, and still 
give the same survival probability (in finite time interval). 

The explanation is that the GF is multi valued, and dif- 
ferent sheets of the function can be used to calculate inte- 
gral ([6]). These ideas can be easily understood for simple 
examples. Consider a site coupled to a semi-infinite lat- 
tice '41 . The system is described by the tight-binding 
Hamiltonian 



^^-^E(l") ('^+11 + 1^ + 1) {n\) 

71=1 

+ 6 Id) {d\-V{\d){l\ + \l){d\), 



(27) 



where \n) is the state localized at the n-th site of the 
lattice. The band (lattice) states are described by the 
Hamiltonian 



Hq = — cos k\k) {k\ 



(28) 



where \k) — ^/2N sin(fcn) \n). Hence we regain 
Hamiltonian ([T|) with 14 = — \/2V^sinfc. After simple 
algebra we obtain (in the upper half plain) 



(29) 



where the square root is defined as having the phase 7r/2 
just above the real axis between —1 and 1, and Aq = 
2V^ . We immediately see that the GF for this model is a 
double valued function, the branch points being +1 and 
— 1. The poles are given by the equation 
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e(l - Aq) ± AoVe^ - 1 2Ao 
1 - 2An 



(30) 



One sheet has real poles for Aq > (e^ -I- l)/2. For Aq — 
(e^ — l)/2 the GF has a second order pole dX uj — (e^ -f 
l)/2e. When Aq increases, this second order pole is split 
into two first order poles, one going right (we assume 
e > 0) and at Aq = 1/2 becoming a pole at the infinity. 
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FIG. 2: Position of poles of the GF for the Hamiltonian ([27] 
for different values of A. The real poles appear on the stan- 
dard sheet, the complex poles on the second sheet. 



For Aq > 1/2 this pole appears for a; < —1. The second 
first order pole, when Aq increases initially approaches 
the point uj = 1, and at a further increase of Aq moves in 
the opposite direction and asymptotically goes to infinity. 

For Ao < (e^ + l)/2 the second sheet has two complex 
poles of the first order. For Aq <C 1 the pole in the 
lower half-plain is situated at e — iAo\/l — and is just 
the FGR pole mentioned above. The poles position is 
presented on Fig. 2. 

In general case the points Eh and Et are the S(ijj) 
branch points. The standard approach consists of making 
the cut along the straight line between the branch points 
and considering only one sheet, thus making the analytic 
continuation into the lower lo half-plain by continuing 
I](w) along the curves which circumvent the right branch 
point clockwise and the and the left branch point anti- 
clockwise On the other hand, we could use a different 
continuation, along the curve passing through the part of 
real axis between Eh and Et- In this case, Eq. ((26|) which 
is valid in the perturbative regime near the point cj = e 
in the upper half-plane is valid in lower half-plain also, 
thus giving a pole in a different sheet of the multivalued 
propagator. 

This pole becomes important if the analytical continu- 
ation of g{uj) into the lower half-plain is done by making 
the cuts from the branch points to infinity and continuing 
the function between the cuts along the curves passing 
through the part of real axis between Eh and £'t, and 
outside as we did it previously. This way to make an an- 
alytic continuation, and hence to calculate the integral 
([6]) is presented on Fig. 3. Due to the exponential de- 
crease of the integrand on the vertical line, in contrast 
to the oscillatory dependence on the real axis, such an- 
alytic continuation is more convenient for the numerical 
calculations for the large t behavior of the non-decay am- 
plitude. 

The second example is defined by the equation 



1^{E) = Ao = const for \E\ < 1. 



(31) 
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FIG. 3: Alternative way to analytically continue the propa- 
gator into the lower half-plane. The arrows show the curves 
of analytical continuation in between the cuts and outside. 
Radius of the big arc and length of the cuts go to infinity. 




FIG. 4: Survival probability as a function of time for Ao — 
.02. 



Thus we get 



Y.{uj) = Ao log 



UJ - 1 



(32) 



The Riemann surface has an infinite number of sheets. 
The standard sheet is obtained by defining log as having 
the phase —it just above the real axis between —1 and 1. 
This sheet always has two real poles, one for w > 1, and 
the other for a; < — 1. 

For the pedagogic purposes let us presents the results 
of numerical calculations for the model considered. The 
time will be measured in units of the FGR time t 



1/r = 27rAo. 



(33) 



For the sake of definiteness we will chose e = —.4. For 
Aq = .02 (see Fig. 4) we observe the FGR regime, say, 
up to i = 9. For A = .1 (see Fig. 5) the FGR regime 
is seen up to t = 3. For Aq = .2 (see Fig. 6) the FGR 



regime is absent. The Rabi oscillations we see already at 
Fig. 5 and still more vividly at Fig. 6. 

In p 




FIG. 5: Survival probability as a function of time for Ao = .1. 
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FIG. 6: Survival probability as a function of time for Ao = .2. 

DEDICATION 

With great pleasure I recall my meetings with I. Vagner 
and our talks both about physics, and life in general. In 
addition to being a very good physicist, he also was a 
very good educator, investing a lot of time and effort in 
his pedagogical activities. So I would like to dedicate this 
my modest contribution, mainly of pedagogical nature, to 
the memory of I. Vagner, who's untimely death I deeply 
regret. 
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